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ABSTRACT 

In a novel approach to studying viscous accretion flows, viscosity has been introduced as a 
perturbative effect, involving a first-order correction in the a- viscosity parameter This method 
reduces the problem of solving a second-order nonlinear differential equation (Navier-Stokes 
equation) to that of an effective first-order equation. Viscosity breaks down the invariance of 
the equilibrium conditions for stationary inflow and outflow solutions, and distinguishes ac- 
cretion from wind. Under a dynamical systems classification, the only feasible critical points 
of this "quasi-viscous" flow are saddle points and spirals. A linearised and radially propagat- 
ing time-dependent perturbation gives rise to secular instability on large spatial scales of the 
disc. Further, on these same length scales, the velocity evolution equation of the quasi-viscous 
flow has been transformed to bear a formal closeness with Schrodinger's equation with a re- 
pulsive potential. Compatible with the transport of angular momentum to the outer regions of 
the disc, a viscosity-limited length scale has been defined for the full spatial extent over which 
the accretion process would be viable. 
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1 INTRODUCTION 

The role of viscosit y in the formation of accretion disc s has, over the years, been recognised to be of paramount importance ( Lvnden-Belll 



196^ : 



Shakura & Sunvaev.,1973 : Pringle 1981: Fr ank et alj|2002h . The standard Keplerian distribution of gaseous matter around a central 



accretor is determined by viscosity, without which there could be no angular momentum transport on outer length scales, and, there fore 



no infall. Viscosity in a Keplerian disc also has a bearing on the time scale of the inward radial drift of matter jprank et al.ll2002l) . So 
viscosity leaves its imprint on accretion processes in more ways than one. While these facts are not a matter of doubt anymore, knowl- 
edge of the exact nature of viscosity still proves elusive. No help i s also forthcomin g from the fact that the observables of an accretion 



disc have been theoretically shown to be independent of viscosity dFrank et al.ll2002r) . To explain the enhanced outward transport of an- 



gular momentum, and t he accompanying inflow rate, it has been variously suggested that turbulence, ordinarily hydrodynamic, or even 



accretion-related studies has been devoted to viscosity from one perspective or the other ('shakura & Sunyaev 


'1973'; 'Liang & Thomson' 1980|; 


Pringle 1981; Matsumoto et al. 1984; Muchotrzeb-Czernv 1986: Abramowicz et al. 1988: Naravan & Yi 


1994; Chakrabarti & Titarchuk 


1995 


Chen et alj 19971: Peitz & AddJ 1997; Frank et al.l2002: Afshordi & Paczvnskj2003 : 


Chakrabarti & Dasl2004l:lBecker & Subramanian 


2005 


Umurhan et al.ll2006l:lDasl 20071: Sharmd 2008: LanzafamelboOSi: Subramanian et al. 


20081) 



In contrast to a viscosity-driven accretion process, a nother model — the sub-Keplerian low angular momentum inviscid flow — has 
by now also become well-established in accretion studies 1 Abramowicz & Zureklll981 ; Fukuelll987 ; Chakrabarti| 1 98S ; Nakavama & Fukue 



I989I: Ichakrabartilll990l : Ixafatos & Ym3ll994 lYang & Kafatoslll 9951 ; IParievlll 9961 ; iMolteni et alJl996l ; lLu et al]|l997l ; lDaj2002l : lDas et al 
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200ilRa\]E003iJ : lBarai etai]|2004lDa3l2004 [Abraham et alJliooj i lPas et ZIboOTi : Ichaudhurv et alJboodlGoswami et alJlioOTi) . This is 
a particularly expedient and simple physical system to investigate, and is considered especially appropriate in describing the the rotating 
flow in the innermost regions of the disc, very close to the event horizon of a black hole. Steady global solutions of inviscid axisymmetric 
accretion on to a black hole have been meticulously studied over the years, and at present there exists an extensive body of literature devoted 
to the subject, with especial emphasis on the transonic nature of solutions, the multitransonic character of the flow, formation of shocks, and 
the stability of global solutions under time-dependent linearised perturbations. 

Having stressed the usefulness of the inviscid model among researchers in accretion astrophysics, it should also be recognised that this 
model has its own limitations. It is easy to understand that while the presence of angular momentum leads to the formation of an accretion 
disc in the first place, a physical mechanism should also be found for the outward transport of angular momentum, especially if its distribution 
is not sub-Keplerian (which, for instance, is the case for strongly-coupled black hole binary systems). This should then make possible the 
inward drift of the accreting matter into the potential well of the accretor. It has already been mentioned that viscosity has been known all 
along t o be a just such a physical means to effect in fall, although the exact presciption for viscosity in an accretion disc is still a matter of much 
debate l Papaloizou & LiJ 1995 : Frank et aLlb002 ). What is well appreciated, however, is that the visco us prescription should be co mpatible 
with an enhanced outward transport of angular momentum. The very well-known a parametrisation of lShakura & SunvaevI ( Il973h is based 
on this principle. 

And so it transpires that on global scales — especially on the very largest scales of a non-sub-Keplerian disc — the inviscid model will 
encounter difficulties in the face of the fact that without an effective outward transport of angular momentum, the accretion process cannot 
be sustained globally. To address this issue, what is being introduced in this paper is t he "quasi-viscous" disc mod el. This model involves 
prescribing a very small first-order viscous correction in the Q-viscosity parameter of IShakura & SunvaevI ( 119730 , about the zeroth-order 
inviscid solution. In doing this, a viscous generalisation of the inviscid flow can be logically extended to capture the important physical 
properties of accretion discs on large length scales, without compromising on the fundamentally simple and elegant features of the inviscid 
model. This is the single most appealing aspect of the quasi-viscous disc model vis-a-vis many other standard models of axisymmetric 
flows which involve viscosity. To dwell on this point further, while Keplerian discs explain infall processes satisfactorily, as far as the 
role of viscosity is concerned, there is the difficulty that the net force driving infall is practically zero (resulting from an exact balance of the 
centrifugal effects against gravity). On the other hand, while sub-Keplerian inviscid discs are free of this difficulty, they do not account for any 
direct outward transport of angular momenturr0 — something that is also necessary to bring about infall. The truth probably lies somewhere 
in between. The quasi-viscous model tries to address that possible area of convergence. While it accounts for an angular momentum transport, 
it also ensures that there is an effective force in the flow to drive the accretion process from an outer boundary to the event horizon of a black 
hole accretor. 

From a most general fl uid dynamical viewpoint, the effect of viscosity is described by a second-order nonlinear differential equation — 
the Navier-Stokes equation jLandau & Lifshitzll987 ). The inviscid limit, on the other hand, is mathematically founded on Euler's equation, 
which is a first-order nonlinear differential equation. The quasi-viscous flow here is based on a perturbative scheme about the inviscid 
conditions, and so the governing equation for this kind of viscous flow can be suitably approximated to a first-order equat ion. In fluid 
dynamics, this is not a particularly unusual mathematical expedient when it comes to accounting for viscosity i Bohr et al]ll993 ). 

The immediate effect of viscosity on stationary flow solutions has been to break down the invariance of the equilibrium conditions 
of inflows and outflows, something that is otherwise preserved well in the inviscid limit. The equilibrium conditions of the quasi-viscous 
flow have been precisely identified, and the nature of the equilibrium points (critical points) has been discerned by devising a first-order 
autonomous dynamical system from the flow equations. In this manner it has been shown that the possible critical points in the phase plot can 
be eit her saddle points, or spi rals, or nodes. In the inviscid limit an earlier study has s hown that only saddle p oints and centre-type points can 
exist I Chaudhurv et al. 2006h . Now centre-type points are a limiting case of spirals ( Jordan & Smith 19991) . Since the quasi-viscous model 
represents a generalisation of the global inviscid flow, but at the same time also implies that viscosity can be tuned to arbitrarily small values, 
a likely scenario that emerges is that the centre-type points (associated with inviscid flows) will become spirals on the inclusion of viscosity 
in the fl ow, however small . This shou ld have various ramifications, especially about connecting multiple transonic solutions through standing 
shocks Jchakrabartilll989l:lDaj2002h. 

An earlier work Jshattachariee & Ray 2007 ) on the quasi-viscous flow, driven by the classical Newtonian potential, has revealed an 
instability — secular instability — when the stationary flow solutions are subjected to small time-dependent perturbations. By analogy, 
exa ctly this kind of insta bility is also seen to develop in Maclaurin spheroids on the introduction of a kinematic viscosity to a first or- 
der I Chandrasekhai 19871) . Similar features in the quasi-viscous flow have also been demonstrated in this paper, but in this instance, under 
the pseudo-Schwarzschild generalisation. 

On large length scales the quasi-viscous flow displays some interesting asymptotic behaviour. Under highly subsonic conditions, the 
pertinent flow equation (Navier Stokes equation) can be transformed mathematically into an equation that resembles Schrodinger's equation 
with a repulsive potential. This has been physically connected to the accumulation of angular momentum on large scales (a property of the 



^ However, in the absence of any viscous t ransport of a ngular momentum, jets launched from accretion discs are supposed to be the only outlet for the intrinsic 
angular momentum of the infalling matter iwiitall200ll) . 
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disc that is very probably also related to the growth of the secular instability on large length scales), and a limiting scale of length has been 
derived from this condition. 

Finally, it should be worth stressing the fact that the entire treatment presented here has been completely analytic, and to the extent that 
this work purports to study pseudo-Schwarzschild flows, it has accounted for the use of any kind of generalised pseudo-Newtonian potential 
to drive the accretion process. This study is the first in a series, in which quasi-viscous accretion around a rotating black hole will also be 
taken up later. This will reveal the influence of black hole spin angular momentum (Kerr parameter) on various equilibrium and stability 
criteria for the flow. 



2 THE EQUATIONS OF THE QUASI- VISCOUS AXISYMMETRIC FLOW AND ITS EQUILIBRIUM CONDITIONS 

For the thin disc, under the condition of hydrostatic equilibrium along the vertical direction iMatsumoto et al.ll 19841 : [Frank et alJbooA . two 
of the relevant flow variables are the drift velocity, v, and the surface density, E. In the thin-disc approximation the latter has been defined by 
vertically integrating the volume density, p, over the disc thickness, H(r). This gives E = pH, and in terms of E, the continuity equation is 
set down as 



aE 1 a ,^ , 

— + -— E«r- 
ot r or 



0. 



(1) 



The axisymmetric accretion flow, driven by the gravitational field of a centrally located black hole, is described in terms of the Newtonian 
geometry of space and time with the help of what is known as a pseudo-Newtonian potential, <j!>(r). This paper will make use of such a general 
expression for the potential, and so the analytical results presented here will hold good under the choice of any pseudo-Newtonian potential. 
Assumption of the hydrostatic equilibrium in the vertical direction will give the condition 

r Ca 

in which the local speed of sound, Cs, and the local Keplerian velocity, vk, are, respectively, defined as = jP/ p and — r4>', with the 
pressure, P, itself being expressed in terms of a polytropic equation of state, P — Kp'' (following this, the speed of sound may also be given 
as Cg — dP/dp). Written explicitly, the disc height can, therefore, be written as 



H 



(2) 



H = {^K) 



and with the use of this result, the continuity equation could then be recast as 



A r„{7+i)/2 



VW d 



p(7+l)/2„j,3/2 



= 0. 



Qtl^ J ' j.3/2 

The general condition for the balance of specific angular momentum in the flow is given by jFrank et alj2002l) . 



27rr \dr J 



dt ^ ' r dr 
where Q. is the local angular velocity of the flow, while the torque is given as 



27rri^Er' 



(3) 



(4) 



(5) 



(6) 



V 9r / ' 

with u being the kinematic viscos ity associated with the flow. With the use of the continuity equation, as equation (O gives it, and going 
by the Shakura&^unvaev 
form jFrank et al.ll2002l : lNM-avan & Yil 19941) 



lll973h prescript ion for the kinematic viscosity, v = aCsH, it would be easy to reduce equation l|5) to the 



d 



V dt 



dr 



pvrH dr 



apHclr-^ ^^Q^. 



with Q,K being defined from vk ~ tJIk. 

Going back to equation a new variable is defined as / 
equation can be closely identified with the constant matter flux rate. In terms of this new variable, equation (O can be modified as 



o(7+l)/2,,^3/2 



A r„{7+l)/21 

dt J 



+ 



r3/2 \dr J 



= 0, 



while equation Q can be rendered as 

d 



ld_ 

V dt 



{r'Q) + 



dr 



/ dr 



dr 



(7) 



/v^, whose steady value, as it is very easy to see from 



(8) 



(9) 



The inviscid disc model is given by the requirement that r^Q, = \, in which A is the constant sp ecific angular momentum. T he quasi- 
viscous disc that is being proposed here will introduce a first-order correction in terms involving a, the Shakura & Sunvaev i 19731) viscosity 
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parameter, about the constant angular momentum solution. Mathematically this will be represented by the prescription of an effective specific 
angular momentum, 



(10) 



with the form of Cl having to be determined from equation {9](, under the stipulation that the dimensionless a-viscosity parameter is much 
smaller than unity. This smallness of the quasi-viscous correction induces only very small changes on the constant angular momentum 
background, and, therefore, neglecting all orders of a higher than the first, a nd ignoring any explicit tim e -variation of the viscous correction 
term, the latter being a standard method adopted also for Keplerian flows i Lightman & Eardlev 1974 ; Shakura & Sunvaev 19761 : Pringle 
198ll:lFranketai]|2002l) . the dependence of f7 on u and p is obtained as 



n = - 



2X-yK 



+ 



With Cl thus defined, it becomes possible under stationary conditions to set down equation l llOb in a modified form as 



Acff(r) = A-2qA ^ 

' VVK 



(11) 



(12) 



Lastly, the equation for radial momentum balance in the flow will also have to be modified under the condition of quasi-viscous 
dissipation. This has to be done according to the scheme outlined in equation jlOj by which, the centrifugal term, \^g{r)/r'^, of the radial 
momentum balance equation, will have to be corrected upto a first order in a. This will finally lead to the result 



dv dv IdP J, , A^ 
at or p or r'^ 



24 {r^^ 



0, 



(13) 



with fl being given by equation dl lb . and P being expressed as a function of p with the help of a polytropic equation of state, as it has been 
mentioned earlier. The steady solution of equation J13t is given as 



dv 1 dP 



A^ 



A^ 



V— + --5- +<?i'(r) - — +4a— ^ 
dr p dr r-^ r-' \ vvk 



0, 



(14) 



whose first integral cannot be obtained analytically because of the Q-dependent term. In the inviscid limit, though, the inte gral is easily 
obtained. This case will be gove rned by conserved conditions, and its solutions have been well-known in accretion literature jChakrabarti 
19891 : IPaslbooi i lPas et al.ll2003h . They will either be open solutions passing through saddle points or closed paths about centre-type points. 
The slightest presence of viscous dissipation, however, will radically alter the nature of solutions seen in the inviscid limit, and it may be 
easily understo od that solutions forming closed paths about cen t re-type points, will, under co nditions of small- viscous correction, be of the 
spiralling kind 1 Liang & Thomsonlll98d : Matsumoto et al.iri984 : Afshordi & Paczvnskill2003l) . This state of affairs is appreciated very easily 
by the analogy of the simple harmonic oscillator. In the undamped state the phase trajectories of the oscillator will, very much like th e 
solutions of the inviscid flow, be either closed paths about centre-type points or open paths through saddle points 1 Jordan & SmithI 19991) . 
With the presence of even very weak damping the closed paths change into spiralling solutions. A more detailed analysis in the regard will 
be carried out in Section[3] 



2.1 The fixed points for polytropic flows 

The pressure, P, is prescribed by an equation of state for the flow 1 Chandrasekhai I939[) . As a general polytropic it is given as P = Kp'^, 
where K is a measure of the entropy in the flow and 7 is the polytropic exponent. In terms of 7, the polytropic index, n, is defined as 
n = (7 — I Chandrasekhai|[l939 ). These definitions are necessary to recast the first integral of equation which is easily obtained for 
stationary conditions. Using the relation between p and Cg, afforded by the polytropic condition, the final expression for the integral could be 
presented as 



2(2n + l)«'r3 



47^2-' ' 



(15) 



where M = ['yK)"m ( Chakrabartill I99(]l ) with m, an integration constant itself, being physically the matter flow rate. 

To obtain the critical points (the equilibrium points) of the flow, it should be necessary to combine both equations l ll4b and U5b . along 
with the polytropic definition of the equation of state, to arrive ultimately at 



2v' 




yv^nj/ 



x' I o2 2 



(16) 



with 0^ — 2(7 + 1)^^. The critical points of the flow will be given by the condition that the entire right hand side of equation ( 116b will 
vanish along with the coefficient of d(u^)/dr. Explicitly written down, following some rearrangement of terms, this will give the two critical 
point conditions as. 
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(17) 



with the subscript "c" labelling critical point values. 

The roots of Tc could be fixed in terms of 7, a, A and M in the r — phase portrait of the stationary flow. In order to do so, the latter 
condition in equation l ll7t can be further modified with the help of the former condition to eliminate v^. This will lead to 



c(/f>'(rc) 



a2 2 
= P 



(18) 



which is a relation that gives rc as a function of 7, a, A and Csc. To eliminate the dependence on Csc, it will be necessary to substitute v in 
equation (HI by using the critical point conditions. This will give 



l/2(n + l) 



(19) 



a result with the obvious implication that the dependence of Tc will finally be given as rc — fi{'y,a,\,A4). One very interesting consequence 
of the presence of viscosity in equation 1 II8I 1 is that the sign of the square root on the right hand side has to be chosen according to whether 
one is studying inflow solutions or outflow solutions. For inflows a negative sign would have to be extracted from the square root, while for 
outflows the chosen sign would have to be positive. This suggests that the invariance of th e coordinates of the cri tical points in the r 



plane would be lost because of dissipation, as opposed to the fully conserved inviscid case Jchaudhurv et ahlbood) . Viscosity, therefore, will 
distinguish accretion solutions from wind solutions. 

The slope of the continuous solutions which could possibly pass through the critical points are to be obtained by applying the L'Hospital 
rule on equation ( I16t at the critical points. This will give a quadratic equation for the slope of stationary solutions at the critical points 
themselves in the r — phase portrait. The resulting expression will read as 



d , 2x 



{Z3 + aZ4) = 0, 



in which the constant coefficients, Zi, Z2, Zz and ^4 are given by 



1 



Z-, 

^3 = — 

7 

and 

^4 = 



7+1 

3-7 



•^"(rc) 



7 ) ^vlrU'(r^) 



— + 20 (rc) + 



3_ _d 

^ ^ sc 1 2 dr y 0' 



8A'd 



respectively. 



7-lMl3_^^^ 



7 + 1 / r, 



(20) 



2.2 The fixed points for isothermal flows 

For an isothermal flow, the appropriate equation of state is given by P = puT / iimn, in which k, is Boltzmann's constant, T is the constant 
temperature, mn is the mass of a hydrogen a tom and /i is the reduced mass, respectively. The definition for H in equation (O will have to 
be modified slightly by setting 7 = 1 ( Afs hordi & Paczvi iski 2003; Ch audhury et al. 2006,). The local speed of sound will also be modified 
to become a global constant of the flow, going by the definition = dP/dp. 

Going back to equation J14b and using the linear dependence between P and p as the appropriate equation of state, will lead to 

v'^+cl^{\np) + ^'{r)~%+Aa%{A\ =0, (21) 
dr dr r'^ r'^ \vvn J 

while the first integral of the stationary continuity condition, given by equation (|4j will give 
As it has been done for polytropic flows, the two foregoing equations can be combined to obtain 
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d , 2^ 27j2 



(2 2\ <-l / 2n 



J,' I -'- 2 



3-r-- 



from which the critical point conditions are easily identified as 



A^ 

rc0'(rc) - -J 

' C 



(23) 



(24) 



Fixing the critical point is much simpler a task for isothermal flows. As it has been done for the polytropic case, Vc has to be eliminated first 
from the second critical point condition in equation l l24t . to obtain 



c0'(rc) 



A^ 



0"(rc) 



<t>'{rc) 



' 0'(rc) 



(25) 



In this expression, the speed of sound, Cg, is globally constant, and so having arrived at the critical point conditions, it should be easy to see 
that Tc and Vc have already been fixed in terms of a global constant of the system. The speed of sound can further be written in terms of 
the temperature of the system as Cs — QT^^^ , where Q — {k/ ^mu)^^^ , and, therefore, it should be entirely possible to give a functional 
dependence for rc, as Tc — f2{a, A, T). The slope of the solutions passing through the critical points in the r — v'^ phase portrait is obtained 
simply by setting 7 = 1 in equation ilOi . 



3 THE CHARACTER OF THE FIXED POINTS: A DYNAMICAL SYSTEMS STUDY 

The s tationary fluid equations describing a viscous flow are in general second-order nonlinear differential equations jLandau & Lifshitz 
1987h . There is as yet no standard prescription for deriving analytic solutions from these equations. Therefore, for any understanding of the 
behaviour of the flow solutions, a numerical integration is in most cases the only recourse. On the other hand, an alternative approach could 
be made to this question, if the governing equations are set up to form a standard first-order dynamical system yprdan & Smith 1999j. The 
mathematical formalism of the stationary quasi- viscous flow is premised on the two first-order differential equations, given by equations ( I14l l 
and J15l >. Of these, the former equation is the result of an approximation to obtain an appropriate first-order differential equation to describe 
a viscous flow. This kind of an approximation is quite usual in general fluid dynamical studies ( Bohr et al. 1993), and short of carrying out 
any numerical integration, this approach allows for gaining physical insight into the behaviour of the flows to a surprising extent. As a first 
step towards this end, for the stationary polytropic flow, as given by equation ( I16t . it should be necessary to parametrise this equation and set 
up a coupled autonomous first-order dynamical system as ( Jordan & Smith 1 99ft) 



2v' 



J,' I /32 2 



dr / 2 

d7 ^ 



(26) 



in which r is a n arbitrary mathematical parameter. With respect to accretion studies in particular, t his kind of parametrisation has been 
repor t ed before (Muchotrzeb-Cze"mvlll986l . iRav & Bhattachariee|[2002l : Ufshordi & Paczvnskill2003l : Ichaudhurv et alJbood : iMandal et aT 



20071 : iGoswami et al.ll2007h . This opens the way to explore the mathematical nature of the critical points much more thoroughly. 

The critical points (which give the equilibrium conditions in the flow) have themselves been fixed in terms of the flow constants. About 



= + 5cf and r = rc + 5r, it becomes 
possible to derive a set of two autonomous first-order linear differential equations in the 5r — ^v^ plane, with itself being expressed in 



these fixed point values, upon using a perturbation prescription of the kind v 



terms of (5r and iv^ from equation ( I15t as 



5c? 
^2 



7 - 



7 + 1 



311 



+ 3 - r-c 



Sr 



(27) 



The resulting coupled set of linear equations in 5r and Sv'^ will be given as 

A5v^ + BSr 



<^ fx 2n 



CSv-" + DSr, 

in which the constant coefficients A, B, C and D are to be read as 

7 - 1 \ v„.2 , 37 - 1 



(28) 



A = 



B 



7 + 1 



7 - 



2y 



, J.I I \ I J.I' I \ I 4>"{^c) 2 T, , 



2 



7+1 
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^«|rj(/)'(r-c) 



7 + 1 



5 rc0"(rc) 
2 2<?!>'(r-c) 



C = 
and 



27 

7 + 1 



7- 



1 \ J. 2 



7+1, 

under the further definition that 



and 



:y = 1 + r, 



<^"'(rc) 



4>"(r,) 



cj>"{r^) (j>'{r,) 



Trying solutions of the type Sv^ ~ exp(57r) and 3r ~ exp(f2r) in equations l l28t . will deliver the eigenvalues 57, which are the growth 



rates of Sv and Sr, as 



{A + D)n + {AD - BC) = 0. 



(29) 



Under a further definition that P — A + D, Q = AD — BC and A = P — 4Q, the solution of the foregoing quadratic equation can be 
written as 



p±Va 



(30) 



Once the position of a critical point, rc, has been ascertained, it is then a straightforward task to find the nature of that critical point 
by using rc in equation ( 130b and all its associated values. Since it has been discussed in Section [2] that rc is a function of a, A and T for 
isothermal flows, and a function of 7, q, A and Ml for polytropic flows, it effectively implies that Q, can, in principle, be rendered as a 
function of the flow parameters for either kind of flow. For isothermal flows, starting from equation ( 123b . a similar expression for the related 
eigenvalues may likewise be derived. The algebra in this case is much simpler and it is easy to show that for isothermal flows the relevant 
results could be derived by simply setting 7 = 1 in equations ( 128b . 

The nature of the po ssible critical points ca n also be predicted from the form of Q in equation ( 130b . If A > 0, then a critical point can be 
either a saddle or a node jjordan&Srniti]|l999h . The precise nature of the critical point will then be dependent on the sign of Q. If Q < 0, 
then the critical point will be a saddle point. Such points are always notoriously u nstable in terms of the sensitivity in generating a solu tion 
through them, after starting from a boundary value far away from the critical point iRav & Bhattacharie el2002l , l2007a|jRov & RavlioOTi) . On 
the other hand, if Q > 0, then the critical point will be a node. Such a point may or may not be stable, depending on the sign of P. If P < 0, 
then the node will be stable. 

A completely different class of critical points will result when A < 0. These points will be like a spiral (a focus). Once again, the 
stability of the spiral will depend on the sign of P. If P < 0, then the spiral will be stable. For inflow solutions in the quasi-viscous disc, the 
form of P (deriving from the sum of A and D) shows that it will always be negative. This is because for inflows the negative root of Vc has 
to be extracted (i.e. Vc < 0) from the square root in the definition of A. Which will obviously mean that if the critical point is either a spiral 
or a node, then it will be stable, with flow solutions in the neighbourhood of the critical point converging towards it. 

The quasi-viscous prescription is based on the requirement that viscosity will only have a small perturbative effect about the conserved 
inviscid flow. In other words, one could tune the viscosity parameter, a, to arbitrarily small (but non-zero) values. In this kind of a situation 
it is much mo re likely than not that A < , and that the stable critical point will be a spiral (nodal points, however, cannot be ruled out 
completely, as lAfshordi & Paczvnskil ( 120030 have shown). Therefore, the most likely picture that emerges as far as the phase portrait of the 
flow is concerned, is that there will be adjacent unstable saddle points and stable spiral points (adj acent points cannot be both stable or 
unstable simultaneously). This argument is in keeping with an earlier study 1 Chaudhurv et al. 20060 on the inviscid disc, where a generic 
conclusion that was drawn about the critical points was that for a conserved pseudo-Schwarzschild axisymmetric flow driven by any potential, 
the only admissible critical points would be either saddle points or centre-type points. For a saddle point, > 0, while for a centre-type 
point, < 0, with being real on both occasions. Noting that a centre-type point is merely a special case (P = 0) of a spiral, introduction 
of viscosity as a small perturbative effect certainly represents a physical generalisation. But with this, what is lost from the phase portrait 
of the flow are homoclinic trajectories connecting a saddle point to itself, or even heteroclinic trajectories connecting two saddle points, 
although one might still argue that heteroclinic paths will exist to connect saddle points with spirals. 

Once the behaviour of all the physically relevant critical points has been understood in this way, a complete qualitative picture of the flow 
solutions passing through these points (if they are saddle points), or in the neighbourhood of these points (if they a re spiral points), can b e 
constructed, along with an impression of the direction that these solutions can have in the phase portrait of the flow 1 Jordan & Smith|[r999l) . 
So what does that imply for multitransonicity, especially about flow solutions which can be generated very far away from the black hole 
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accre tor to reach its event horizon eventually? Many earlier studies i Chakrabarti 19891 199d : Daj 2002 ; Das et al. 2003 ; Chaudhurv et al 



2006h have taken up this question in great detail, and it has been shown that for certain parameter-space values pertaining to the inviscid disc , 



three critical points can result. These are located in such a manner that a centre-type point is flanked by two saddle points i Chakrabartil 1 99(]|) 
through which transonic solutions pass. For very small values of viscosity, it is no w conceivable that the cen tre-type point in the middle will 
become a stable spiral. This is in fact very much in keeping with the conclus ion of Liang^^Thomso that the number of independent 

transonic solutions cannot exceed one plus the number of spiral singularities. Liang & TTjomsorJIigS ff) have further suggested that in realistic 
physical situations, models with spiral singularities are unstable and the critical transonic solution wheneve r it exists is unique in rel evant 
situations. While the instability of the quasi-viscous disc will indeed be verified in Section|4l the argument of Liang & ThomsonI ( 1980l) also 
conceivably has a stron g bearing on another feature that is very intimately connected to multitransonicity in accretion flow — shocks, with 
or without dissip ation aChakrabartill 19891 Il99d lDasl2002l : lOas et aklbooa Ichakrabarti & Dasll2004l : lDasl2007l: Ipukumura & Kazanasll 20071 : 
LanzafamdboOSh . 



4 TRAVELLING- WAVE PERTURBATIVE ANALYSIS : SECULAR INSTABILITY 



Many earlier works have tak en up the q uestion of the stability of viscous thin disc accretion i Lightman & Eardlevlll974 ; Shakura & Sunyaev 

1976 : lLivio & Shavivll977l :lKato 1978: Umurhan & Shavivl20o3^ . Regarding quasi-viscous accretion discs in particular, an earlier study dBhattachariee & Ravi 



2007D has shown that stationary flow solutions driven by the simple Newtonian potential suffer from an instability under time-dependent 



perturbations — both as a standing wave and as a high-frequency travelling wave. It need not always be true that standing and travel- 
ling waves will simultaneously exhibit the same qualitative properties as far as stability is concerned. Many i nstances in fluid dynamics 
bear this out. In the case of binary fluids, standing waves indicate instability, as opposed to travelling waves ( Cross & HohenbergI 1993 ; 



■Bhattacharva & Bhattach arjee 200 5J, while the who l e physical picture i s quite the opposite for the fluid dynamical problem of the hydraulic 
jump teohr et aljl993l : Rav & Bh attacharieelboOTbl : Singha et al.ll2005 ). Contrary to all this, the axi symmetric stationary quasi- viscous flow 
is greatly disturbed both by a standing wave and by a travelling wave i Bhattacharjee & Rav 2007h . This provides convincing evidence of 
its unstable charact er, and it is very muc h in consonance with similar conclusions drawn from some earlier studies. For high-frequency 
radial perturbations Chen & TaamI J 19931) have found that inertial-acousti c modes are local ly unstable, with a greater degree of growth for 

1988h have revealed a growth in the amplitude of a 
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the outward travelling modes than the inward ones. On the other hand, 
non-propagating perturbation at the critical point, which, however, becomes stable in the inviscid regime. 

This kind of instability — one that manifests itself only if s ome dissipative mech anism (viscous dissipation in the case of the quasi- 
viscous rotational flow) is operative — is called secular instability i ChandrasekhaJ 1987 ). It should be very much instructive here to furnish a 
parallel instance of the de stabilising influence o f viscous dissipation in a system undergoing rot ation: that of the effect of viscous dissipation 
in a Maclaurin spheroid jChandrasekhaii 1 1 987h . In studying ellipsoidal figures of equilibrium. IChandrasekhad ( Il987h has discussed that a 
secular instability develops in a Maclaurin spheroid, when the stresses derive from an ordinary viscosity which is defined in terms of a 
coefficient of kinematic viscosity (as the a parametrisation is for an accretion disc), and when the effects arising from viscous dissipation 
are considered as small perturbations on the inviscid flow, to be taken into account in the first order only. It is exactly in this spirit that 
the "quasi-viscous" approximation has been prescribed for the thin accretion disc, although, unlike a Maclaurin spheroid, an astrophysical 
accretion disc is an open system. 

Curiously enough, the geometry of the fluid flow also seems to be having a bearing on the issue of stability. The same kind of study, 
as has been done here with v iscosity in a rotational flow, had also been done earlier for a viscous spherically symmetric accreting system. 
In that treatment (IRavll2003bl) viscosity was found to have a stabilising influence on the system, causing a viscosity-dependent decay in the 
amplitude of a linearised standing-wave perturbation. This is quite in keeping with the understanding that the respective roles of viscosity are 
at variance with each other in the two distinctly separate cases of spherically symmetric flows and disc flows. While viscosity contributes to 
the resistance against infall in the former case, in the latter it aids the infall process. 

An important aspect of the time-dependent perturbative analysis presented here is that the flow has been modelled to be driven by a gen- 
eral gravitational potential, (j) (as opposed to the choice of any particular kind of mathematical f orm for 4>). Therefore, none of conclusions re- 
garding secular instability w ill be qualified upon using any of the pseudo-Newtonian potentials ( Paczvnski & Wiitall980l : Nowak & Wagonei 



1991 ; Artemova et al. 199a), which are regularly invoked in accretion-related literature to describe rotational flows on to a Schwarzschild 



black hole. This shall be especially true for the flow on large scales, where all pseudo-Newtonian potentials converge to the Newtonian limit. 
The following treatment will bear this out. 

To proceed with the perturbative analysis, a time-dependent perturbation is introduced about the stationary solutions of the flow vari- 
ables, V and p, according to the scheme, v{r,t) = vo{r) + v'{r,t), p(r,t) — po{r) + p'{r,t) and f{r,t) — fo{r) + f'{r,t), in all of 
which, the subscript "0" implies stationary values. In particular, the stationary solution /o is a constant, as both equations (|4j and ^ indicate. 
This constant, as it is immediately evident from a look at equation l|4j, is very much connected to the matter flow rate, and, therefore, the 
perturbation /' is to be seen as a disturba nce on the steady, cons tant background accretion ra te. For spherically symmetric flows, this Eulerian 
perturbation scheme has been applied by Petterson et al. 1 198C ) and Theuns & David 1 1992 ), while for inviscid axisymmetric flows, the same 
method has been used equally effectively by Ravl i 2003a ) and Chaudhurv et al. 1 2006h ^ 

The definition of / will lead to a linearised dependence among /', v' and p' as 
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Jo V 2 / po vo 

while from equation Q, an exclusive dependence of p' on /' will be obtained as 
0, 



dt '^^ fo \ dr 

with — 2(7 + 1)^^. Combining equations ( 131b and i32i will render the velocity fluctuations as 

dt fo \dt dr 

which, upon a further partial differentiation with respect to time, will give 
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From equation l ll3t the linearised fluctuating part could be extracted as 
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dt dr 
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in which a — Cso/(iio«k) and Cso is the local speed of sound in the steady state. Differentiating equation l l35t partially with respect to t, and 
making use of equations J32t . J33t and ([34} to substitute for all the first and second-order derivatives of v' and p', will deliver the result 
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^ 0, (36) 



entirely in terms of /'. This is the equation of motion for a perturbation imposed on the constant mass flux rate, fo, and it shall be important 
to note here that the choice of a driving potential, Newtonian or pseudo-Newtonian, has no explicit bearing on the form of the equation. 

With a linearised equation of motion for the perturbation having been derived, a solution of the form /'(r, t) = g^{r) exp(— ici;t) is 
applied to it. From equation ([36), this will give 



2 d , , 1 d 
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0. (37) 



The perturbation is now made to behave in the manner of a radially travelling high-frequency wave, whose wavelength is suitably con- 
strained to be small, i.e. it is to be smaller than any characteristic length scale in the system. Effectively this invokes the WKB approximation, 
and a perturb ative treatment of this nature has be en carried out before on spherically symmetric flows ( Petterson et al.|[l980,) and on axisym- 
metric flows JiavlboOSijJchaudhurv etallbooei) . In both these cases the radius of the accretor was chosen as the characteristic length scale 
in question, and the wavelength of the perturbation was required to be much smaller than this length scale. In this study of an axisymmetric 
flow driven by the gravity of a black hole, the radius of the event horizon could be a choice for such a length scale. As a result, the frequency, 
uj, of the waves should be large. 

An algebraic rearrangement of terms in equation ( I37t will lead to an integro-differential equation of the form 



T 



(dr) 



dr = 0, 



(38) 



with its coefficients being given by 

V = vl- (fclo, 
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At this stage, bearing in mind the constraint that uj is large, the spatial part of the perturbation, guj{r), is prescribed as gi^{r) = exp(s) 
with the function s itself being represented as a power series of the form 



= y Ml) 



(39) 



The integral term in equation ( I38t . can, through some suitable algebraic substitutions, be recast as 
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with S being given by another power series as 



m — 1 



(40) 



Following this, all the terms in equation J38b can be expanded with the help of the power series for gui{r). Under the assumption (whose 
self-consistency will be justified soon) that to a leading order 



^ do- ^ do- 

ds dr V dr / 



the three successive highest-order terms (in a decreasing order) involving lu will be obtained as co^, co and a;". The coefficients of each of 
these terms are to be collected first and then individually summed up. This is to be followed by setting each of these sums separately to zero, 
which will yield for uj^, uj and uj'^, respectively, the conditions 
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Out of these, the first two, i.e. equations and l l42t . will deliver the solutions 



va ± /3cso 



dr 
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■ dr, 



(43) 



(44) 
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respectively. 

The two foregoing expressions give the leading terms in the power series of guj{r). While dwelling on this matter, it will also be 
necessary to show that all successive terms of s(r) will self-consistently follow the condition iij~"\k„{r)\ ^ a;~*-"''''^^|A:„+i(r)|, i.e. the 
power series given by g^, (r) will converge very quickly with increasing n. In the inviscid limit, this requirement can be shown to be very 
much true, considering the behaviour of the first three terms in fc„(r) from equations i44\ . J45b and ( 143 1 . These terms can be shown to go 



asymptotically as fc_ 



r, fco 



In r and ki 



. given the condition that no 



5/2 



on large length scales, while Cso approaches its 



constant ambient value. With the inclusion of viscosity as a physical effect, it can be seen from equations l l44t and ([45}, respectively, that while 
k-i remains unaffected, ko acquires an a-dependent term that goes asymptotically as r. This in itself is an indication of the extent to which 
viscosity might alter the inviscid conditions. However, since a has been chosen to be very much less than unity, and since the wavelength of 
the travelling waves is also very small, the self-consistency requirement still holds. Therefore, as far as gaining a qualitative understanding 
of the effect of viscosity is concerned, it should be quite sufficient to truncate the power series expansion of s(r), after considering the two 
leading terms only, and with the help of these two, an expression for the perturbation may then be set down as 
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(46) 



which should be seen as a linear superposition of two waves with arbitrary constants A+ and Both of these two waves move with a 
velocity /3cso relative to the fluid, one against the bulk flow and the other along with it, while the bulk flow itself has a velocity vo. It should 
be immediately evident that all questions pertaining to the growth or decay in the amplitude of the perturbation will be crucially decided 
by the real terms delivered from fco. The viscosity-dependent term is especially crucial in this regard. For the choice of the lower sign in 
the real part of /' in equation ( I46t . i.e. for the outgoing mode of the travelling wave solution, it can be seen that the presence of viscosity 
causes the amplitude of the perturbation to diverge exponentially on large length scales, where Cso — Cs{oo) and vo ~ r"^''^, with —vq 
being positive for inflows. The inwardly travelling mode also displays similar behaviour, albeit to a quantitatively lesser degree. It is an easy 
exercis e to see that stability in the system would be restored for the limit of a = 0, and this particular issue has been discussed bv iRav 
( 2003ah and lchaudhurv et al.l J2006h . The exponenti al growth behaviour o f the amplitude of the perturbation, therefore, is exclusively linked 
to the presence of viscosity. Going back to a work of Chen & Taaml 1 1993 ), it can be seen that the inertial-acoustic modes of short wavelength 
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radial perturbations are locally unstable throughout the disc, with the outward travelling modes growing faster than the inward travelling 
modes in most regions of the disc, all of which is very much in keeping with what equation l l46t indicates here. 
With the help of equation i32t it should be easy to express the density fluctuations in terms of /' as 



and likewise, the velocity fluctuations may be set down from equation iiW as 

^/ ^ A _ «o ds\ f_ 
VQ \ iLu Ar) fo ' 

In a unit volume of the fluid, the kinetic energy content is 

Skill = 7^ {po + p) (yo +v'Y , 
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(48) 



(49) 



while the potential energy per unit volume of the fluid is the sum of the gravitational energy, the rotational energy and the internal energy. 
For a quasi-viscous disc, to a first order in a, this sum is given by 
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where e is the internal energy per unit mass i Landau & LifshitJl987 ). In equation ( 150b the effective angular momentum for the quasi-viscous 
disc will have to be set up as a first-order correction about the inviscid conditions. Following this, a time-dependent perturbation has to be 
imposed about the stationary values of v and p. All first-order terms involving time-dependence in equations i49\ and i50\ will vanish on 
time-averaging. In this situation the leading contribution to the total energy in the perturbation comes from the second-order terms, which 
are all summed as 
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(51) 



In the preceding expression all terms involving p and v' can be written in terms of /' with the help of equations ( 1471 ) and l l48t . in both of 
which, to a leading order, s ~ Lok-i. This is to be followed by a time-averaging over /' , which will contribute a factor of 1/2. The total 
energy flux in the perturbation is obtained by multiplying fpcrt by the propagation velocity {vq ± /3cso) and then by integrating over the area 
of the cylindrical face of the accretion disc, which is 2nrH. Under the thin-disc approximation, H <^ r, this will make it possible to derive 
an estimate for the energy flux as 
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in which M is the Mach number, defined as M = uq/cso, while 
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When M — > on large length scales, ( converges to a finite value. However, on these same length scales, what will not converge are the two 
terms involving a in equation ( 1521 ). Under the asymptotic conditions on vq and Cgo, discussed earlier, one term will diverge exponentially, 
while another will have a power-law growth behaviour of r^. The quasi-viscous disc will, therefore, be unstable on large scales under the 
passage of a linearised radially propagating high-frequency perturbation. It is easy to check that under inviscid conditions, with a = 0, 
and for an adiabatic perturbation with p = 1, the disc will immediately revert to stable behaviour, a feature that has been mathematically 
demonstrated in some earlier papers Ray 2003a : Chaudhurv et al. 20061 ). That instabilities should develop exponentially on large length 
scales, because of the merest presence of viscosity (to a first order in a, which itself is much less than unity) has disturbing implications, 
since all physically meaningful inflow solutions will have to pass through these length scales, connecting the outer boundary of the flow with 
the surface of the accretor (or the event horizon, if the accretor is a black hole). 



5 ASYMPTOTIC BEHAVIOUR OF THE QUASI- VISCOUS DISC 



A previous study dRav & Bhattacharieg|2007ah has dwelt on how it should be possible to select the stationary transonic solution of the 
inviscid axisymmetric flow through a non-perturbative time-dependent criterion (such as the minimisation of the total specific energy of the 
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flow). This kind of insight into the long-time behaviour of the inviscid flow on large length scales of the thin accretion disc allows for making 
a similar foray into the quasi-viscous disc. However, in this instance, a straightforward application of the methods employed for the inviscid 
disc would not be possible on two counts. First, the quasi-viscous disc being a dissipative system (i.e. with energy being allowed to be 
drained away from this system), there should be no occasion to look for the selection of a particular solution, and a selection criterion thereof, 
on the basis of energy minimisation. Secondly, the fact that the quasi-viscous disc is unstable on large length scales, is reason enough to 
believe that no solution — transonic or otherwise — might be free of time-dependence. Therefore, a long-time evolution of the quasi-viscous 
disc towards a stationary end is not something that might be hoped for. For all that, however, this kind of a disc system does exhibit some 
interesting asymptotic features on large length scales. 

On such length scales of an accretion disc, all pseudo-Schwarzschild flows converge to the Newtonian limit, i.e. <j){r) ~ —GM/r. 
Meanwhile the stationary solution of equation Q can be expressed on the same length scales simply as, 

P^,'^+''/'^or-'/^ ~ -mx/GM, (53) 

with m being the conserved matter inflow rate. The negative sign arises because for inflows, vq goes with a negative sign. From equation l l53t . 
with po approaching a constant ambient value on large length scales, the drift velocity, vq, can consequently be seen to go asymptotically 
as r~^^'^ . Bearing in mind that for inflows, vo < 0, the asymptotic dependence of the effective angular momentum can be shown from 
equation l ll2t to be 



Aeff (r) ~ A 



1 + 2a 



id' 



(54) 



where ri is a scale of length, which, to an order-of-magnitude, is given by rf ~ GMm[Cs (cxj)p(oo)]~^. This asymptotic behaviour is 
entirely to be expected, because the physical role of viscosity is to transport angular momentum to large length scales of the accretion disc. 

Nevertheless, the distribution of matter in a viscous disc takes place on a time scale determined by viscosity, and therefore a study of the 
time-dependent properties in a viscous disc is one of the few means of acquiring some im pression about the role of viscosity, especially since 



the observables in a steady disc are largely independent of viscosity i Frank et al]|2002 ). With that objective in mind, it will be worthwhile 



first to try to understand the structure of the governing time-dependent differential equation for the flow on large length scales. To do so, 
it shall be necessary to invoke the approximation that very far from the accretor in the outer regions of the flow, on highly subsonic scales 
of velocity, the density variations are negligibly small compared to the time evolution of the velocity field. The evolution will consequently 
follow the general Navier-Stokes equation in the limit of djVj — 0, which can be set down as 

dtVi + VjdjVi = vdjdjVi — diV, (55) 

where the potential function V = V(r, t) = nc^ + 4>{r) + /2r^ , with n being the usual polytro pic index. Quite evidently, equation 
is a nonlinear differential equation, but with the help of the Hopf-Cole transformation jRegevlbood) , 



v^ = - jftC, (56) 
it can be reduced to a linear form in the variable ^, going as 

2udti = 2//2fta,C + Vi. (57) 

The potential function V will in general be modified by the addition of an integration constant, which might physically be identified from 
the ambient conditions of the fluid. This, however, will require a knowledge of the boundary conditions for the scalar function ^, something 
whose inherent difficulties would be appreciated soon. In the outer regions of the disc, however, the speed of sound, which is a scalar function 
of the density, would asymptotically assume a constant value. 

An interesting aspect of equation ( 1571 ) is that for the one-to-one correspondence of 2u with iTi, and of m with 1, there is an exact 
equivalence between this equation and Schrodinger's equation, 

iUdtip = -—d^di^ + Vip. (58) 
2m 

For the steady limit of the potential function V (a requirement that is satisfied on large length scales), it is easy to carry out a separation of 
variables in equation ( 1571 ), and the resulting stationary eigenvalue equation in ^ may then be expressed in a Hamiltonian form as 

- 2i.V^C + Vi^ E^, (59) 

in which V — —V. A comparison between equation ( |59b and the stationary Schrodinger equation ought to be instructive and insightful here. 
To have any notion of how ^ evolves in time, one would have to determine the eigenvalues given by E, and it is here that a great stumbling 
block is encountered. To solve the second-order differential equation given by equation ( |59t , two boundary conditions are imperative, and 
they in turn would characterise the eigenvalues. The outer boundary condition on ^ is relatively easy to prescribe. For r — > oo and v — > 0, 
the scalar function ^ will asymptotically approach a constant value, and there is much similarity in this with the asymptotic behaviour of 
another scalar function, the density p. Knowing the precise inner boundary condition on 5, however, is a most difficult problem. First of all, 
it depends on the nature of the accretor itself. While a black hole will have an event horizon, a neutron star or an ordinary star will have 
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a physical surface. This fact alone has much influence over the inner bound ary condition. Apart from this, realistically speaking, various 
astrophysical processes near the surface of the accretor will affect the flow ( Petterson et al. 198Ci) . In any case, equation ( |59t is valid for 
large scales only. Yet, short of actually having to solve equation l |59K it should still be possible to derive some information on how viscosity 
affects the flow in the outer regions of the disc. It is important to see that V in equation ( I59t assumes the properties of a repulsive potential. 
More than the pressure effects, this repulsive nature is reflective of the cumulative transfer of angular momentum to large length scales of an 
accretion disc. Referring to equation l l54b . one can see that in the outer regions of the disc, the effective specific force, (p, is given to a first 
order in a by 



GM 



73 +''"73"' 



(60) 



from which it is evident that on scales of r ~ ri, the transport of angular momentum will give rise to an asymptotic constant non-zero force 
opposed to gravity. For the viscous disc this gives rise to a repulsive effect on large length scales, in opposition to gravitational attraction. 
From this argument one may go a step beyond and conjecture that ria~^^^ defines a limiting length scale for accretion, that the outward 
transport of angular momentum imposes. This is also the length scale on which secular instability is most conspicuous. In units of the 
Schwarzschild radius of the black hole, 2GM/(? , this length scale is seen to be 

1/3 

(61) 
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All of this is quite compatible with how viscosity redistributes an annulus of matter in a Keplerian flow around an accretor; the inner region 
drifting in because of dissipation, and consequently, through the conservation of angular momentum and its outward transp ort, making it 
necessary for the outer regions of the matter distribution to spread even further outwards I Pringle 1981 ; Frank et al. 2002h . This state of 
affairs is qualitatively not altered in anyway for the quasi-viscous flow, except for the fact that with viscosity being very weak here, the 
outward transport of angular momentum can perceptibly cause an outward drift of matter only on very large scales. It is obvious that once 
Q = 0, i.e. for the inviscid limit, as equation | |6U shows, this scale would be shifted to infinity. 



6 CONCLUDING REMARKS 



One very important physical role of viscosity in an accretion disc is that it governs the distribution of matter in the disc. The manner in which 
viscosity redistributes an annulus of matter in a Keplerian flow around an accretor is very well known, with the inner region of this disc 
system drifting in because of dissipation, and consequently making it necessary for the ou ter regions of the matter distribu tion to spread out 
even farther, because of the conservation of angular momentum and its outward transport ( |Pringlelll98ll : iFrank et al I2OO2I) . 

Viscosity, however, also gives rise to secular instability in the quasi-viscous disc. This casts much doubt on the long-term viability of 
the accretion flow, and its temporal evolution towards a stationary state. It may rightly be argued that the instability that develops on the 
large subsonic scales of a quasi- viscous disc is intimately connected with the cumulative transfer of angular momentum on these very length 
scales. The accumulation of angular momentum in this region may create an abrupt centrifugal barrier against any further smooth inflow of 
matter. However, this adverse effect could disappear if there should be some other means of transporting angular momentum from the inner 
regions of the disc. Astrophysical jets could readily afford such a means Iwiitai2001\ insofar as jets actually cause a physical drift of angular 
momentum vertically away from the plane of the disc, instead of along its radial length. This will be all the more true if this off-the-plane 
angular momentum drift happens on length scales that are much smaller than the scale indicate d by e quation ( 161b . 

Stability could be restored through many other means. A recent work bv lMach & Maled ( 1200 8) has shown numerically how it should 
be possible to have stable steady accretion solutions for tr ansonic flows of a self-gravi tating gas. This stability argument holds true even 
for perturbations in the nonlinear regime. Another work by Nagakura & Yamad3 1 200 sh has established the stability of an accretion shock 
(connecting transonic solutions passing through two distinct saddle points) under axisymmetric perturbations. 
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